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We study conductance fluctuations in disordered quantum wires with unitary symmetry 
focusing on the case in which the number of conducting channels in one propagating direction 
is not equal to that in the opposite direction. We consider disordered wires with N -\- m left- 
moving channels and A'' right-moving channels. In this case, m left-moving channels become 
perfectly conducting, and the dimensionless conductance g for the left-moving channels behaves 
as (? ^ m in the long- wire limit. We obtain the variance of g in the diffusive regime by using 
the Dorokhov-Mello-Pereyra-Kumar equation for transmission eigenvalues. It is shown that the 
universality of conductance fluctuations breaks down for m 7^ unless A*' is very large. 
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1. Introduction 

I The discovery of a perfectly conducting channel in dis- 
' ordered wires^"^'' provides a counterexample to the con- 
[ jecture that an ordinary quasi-one-dimensional quantum 
' system with disorder exhibits Anderson localization (i.e., 
' conductance decays exponentially with increasing system 

length L and eventually vanishes in the hmit of L ^ 00) . 

We have shown that perfectly conducting channels can 

• be stabilized in two standard universality classes. One 
I is the symplectic universality class with an odd number 

■ of conducting channels, and the other is the unitary 
I universality class with the imbalance between the num- 
1 bers of conducting channels in two propagating direc- 
I tions.^"^^ The symplectic class consists of systems hav- 

■ ing time-reversal symmetry without spin-rotation invari- 
I ance, while the unitary class is characterized by the ab- 

■ sence of time- reversal symmetry. 

Much attention has recently been paid to electron 
transport in the above-mentioned two universality classes 
] which do not exhibit Anderson localization. For the sym- 

• plectic class, one perfectly conducting channel is stabi- 

■ lized in the odd-channel case, while such a special channel 

• does not exist in the ordinary even-channel case. We have 
' studied in details how the dimensionless conductance g 

behaves as a function of LJ'^ It is shown that the be- 
havior of g in the odd-channel case is very different from 
that in the even-channel case when L is much longer than 
the conductance decay length ^. The dimensionless con- 
ductance in the even-channel case decays as g — *■ with 
increasing L, while 17 ^ 1 in the odd-channel case due 
to the presence of one perfectly conducting channel. Fur- 
thermore, the decay of g with increasing L is much faster 
in the odd-channel case than in the even-channel case. 
However, such a notable even-odd difference does not 
appear in the diffusive regime of L <C ^.^'^ For the uni- 
tary class, the number of perfectly conducting channels 
depends on the channel-number imbalance between two 
propagating directions. If m perfectly conducting chan- 
nels are present, the dimensionless conductance behaves 



as (7 — > m with increasing L. The present authors have 
shown that the notable m-dependence of g appears in the 
long-L regime, by using a scaling approach^^ and a super- 
symmetry approach. However, the behavior of g in 
the diffusive regime has not been clarified. It is interest- 
ing to study whether the perfectly conducting channels 
affect the behavior of g in the diffusive regime. 

In this paper, we focus on the unitary universality 
class with the channel-number imbalance and consider 
the conductance in the diffusive regime. We here present 
the basic framework to describe the electron transport 
in systems with the channel-number imbalance on the 
basis of the scaling approach. Let us consider disor- 
dered wires with N right-moving channels and N -\- m 
left-moving channels. In this case, m left-moving chan- 
nels become perfectly conducting, and the dimensionless 
conductances g and g' for the left-moving and right- 
moving channels, respectively, differ from each other. 
Each dimensionless conductance is determined by a cor- 
responding set of transmission eigenvalues. If the set of 
the transmission eigenvalues for the right-moving chan- 
nels is {Ti, T2, . . . , Tn}, that for the left-moving channels 
is expressed as {Ti, T2, . . . , Tn, 1, . . . , 1}, where we have 
identified the iV -f 1 to -I- mth channels as the per- 
fectly conducting ones. The dimensionless conductance 
g is given hy g = Yfl^^ Ta = m + YlH^i^a, while 
g' — Y^^=i '^a- We observe that g = g' + m. It should 
be noted that the mean free paths I and I' for the left- 
moving and right-moving channels, respectively, are not 
equal due to the channel-number imbalance. Indeed, they 
satisfy I' — {N/{N -f m))l. The statistical behavior of g, 
as well as g' , is described by the probability distribu- 
tion for {Ti,T2, . . . ,Tjv}. We define A, = (1 - Ta)/T, 
and introduce the probability distribution P{{Xa}; s) for 
{Ai, A2, . . . , Xn}, where s is the normalized system length 
defined by s = L/l. The Fokker-Planck equation for 
P({Aa}; s), which is called the Dorokhov-McUo-Pereyra- 
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Kumar (DMPK) equation/^' is expressed as^^ 
dP{{Xa};s) 



function F{{Xa}) is defined by 



ds 

1 ^ 
N ^ 



d 



Aa(l + A„)J 



_d_ fPiiXahs] 
dXa 



with 



J 



N N-1 N 

u^Txu n 1^-^ 

c=l b=l a=b+l 



J 



(1) 



(2) 



The factor H^i A™ in J represents the repulsion aris- 
ing from the m-fold degenerate perfectly conducting 
eigenvalue. This eigenvalue repulsion suppresses the 
non-perfectly conducting eigenvalues {Ti, T2, . . . , TV}- It 
should be emphasized that all the influences of the per- 
fectly conducting channels are described by this factor. 

The purpose of this paper is to study how the per- 
fectly conducting channels affect the dimensionless con- 
ductance g in the diffusive regime. Particularly, we focus 
on the variance, var{(/} = (g^) — (5)^, which character- 
izes the magnitude of conductance fluctuations. We know 
that in the ordinary case of m = 0, the variance takes 
the universal value 1/15 irrespective of the normalized 
system length s = L/l. This is called universal conduc- 
tance fluctuations. Does this universality hold even in the 
presence of the perfectly conducting channels? To answer 
this question, we obtain vax{g} by using two approaches 
based on the DMPK equation. First, we analytically cal- 
culate the variance as a function of s by using an A''"^ ex- 
pansion approach. This is applicable to the case in which 
N » max{m, 1}. We show that the variance is given by 
varj^} = 1/15 irrespective of m in the diffusive regime of 
N ^ s ^ 1. This indicates that in the large- TV limit, the 
universality of conductance fluctuations hold even in the 
case of m 7^ 0. To study the case in which the ratio m/N 
is not very small, we numerically calculate varj^"} by us- 
ing a classical Monte Carlo approach based on an approx- 
imate probability distribution for transmission eigenval- 
ues. This is applicable to the case of an arbitrary m as 
long as s/4:N <gc 1. We treat the cases of N = 6 and 20 
with TO = 0, 1, 2, 3. We show that in the case of to = 0, 
the variance approximately takes a constant value nearly 
equal to 1/15 for > s > 1. However, deviation from 
this universal behavior arises when to ^ 0. We show that 
for TO ^ 0, the variance does not take a constant value, 
but decreases with increasing s, and the corresponding 
deviation becomes more pronounced with increasing to. 
We also show that the deviation in the case of TV = 6 is 
more noteworthy than that in the case oi N = 20. These 
results indicate that the universality of conductance fluc- 
tuations breaks down for m 7^ unless N is very large. 

In the next section, we analytically obtain varj^} by 
using the expansion approach. In §3, we numeri- 

cally calculate varj^} by using the classical Monte Carlo 
approach. The Monte Carlo results are compared with 
those obtained in §2. Section 4 is devoted to summary. 

2. Expansion Approach 

We consider var{g'} = (ff^) — (5)^ in the diffusive regime 
oi N ^ s ^ 1. The ensemble average of an arbitrary 



(F)^ dAi---dAwF({Aa})P({A4;s). (3) 
Jo 

Let F be 

N N 

r-Er.^E^. (4) 

a— 1 a—1 

We can express varj^i} in terms of (F^) with p = 1 and 
2. To obtain (F^*), we derive the evolution equation for 
it on the basis of the DMPK equation. Combining the 
DMPK equation with eq. (3), we obtain 

d{F) 



N- 



{|:7^{-<--)^l£}> 

I ^ ( rP'F r)F 

(g{A.(I+A.)^ + (I+2AJ^ 

/ N 

+ A„(1 + A„) T 

V b=i ■^'^ 



_^m\dF' 

Xb Xa I dXa 



(6#o) 

Replacing F by F^, we obtain 



(5) 



iV- 



ds 



= -mp{TP 



+ pip - l)(FP-2 (F2 - F3)) -p{TP+'), (6) 



where 



(7) 



Note that cq. (6) is not closed since its right-hand side 
contains (PP-^Fz) and (PP-^Fg). Even if eq. (6) is com- 
bined with differential equations for (r^'r2) and (F^Ts), 
they do not form a closed set of equations and we can- 
not obtain (P^) in a simple manner. To overcome this 
diSiculty, we adapt the expansion approach pre- 

sented by Mello and Stone.-'* This approach is adaptable 
if » max{TO, 1}. We expand (P^) in a power series of 
N~'^, and obtain it up to order of N~'^~^p. To do so, we 
must supplement eq. (6) by the following equations, 

TV^^^^ = -mip + 2)(rf Fa) + M'^'"' i^s - ^i)) 

+ p(p-l)(FP-2 (r^-PsPg)) 

+ 2(Ff+2)-(p + 4)(FP+iF2), (8) 

^^^^7^ = -m(p + 3)(Ff F3) + 6p{TP-' (F4 - Fs)) 
+ p{p -i){rp-' {r.Ts-Tl)) 
+ 6(Ff+iF2)-(p + 6)(FP+iF3) 

-3(FJTi), 



(9) 



N 



ds 



TO(p + 4)(FPr2)+8(FP(F4-F5)) 



J. Phys. Soc. Jpn. 



Full Paper 



Yositcike Takane and Katsunori Wakabayashi 



+ 8p(rf-i(r2r3-r2r4)) 
+p{p -i){rp-^ (r^-r^ra)) 

+ A{TP+^T2) - (p + 8){TP+^Tl). (10) 

These equations can be derived from eq. (5). We employ 
the following expansions, 

{TP) = NPfp^ois) + A^^-Vp,i(s) + NP-^fp,2{s) + ■■■, 

(11) 

(rfr2) = 7VP+Vp+i,o(s) + ---, (12) 

{TPT3) = NP+'k^+Us) + -.- , (13) 

{rprl) = NP+%+2,o{s) + --- (14) 

with 

/p,n(0) = Jp,n(0) = kp,n{s) = lp,n{0) = <5„,0- (15) 

Substituting these expansions into eqs. (6) and (8)-(10), 
and equating the coefficients of the various powers of N, 
we obtain a set of closed differential equations for /p,o(s), 
fp,i{s), fpAs), ip,o(s), fcp,o(s) and lpfl{s). We solve the 
resulting equations under the initial conditions given in 
eq. (15). These procedures are outlined in Appendix. We 
finally obtain 



+ 



+ 



NP 

JTTsjp ^ 2(1 + s)P+i 



[s' + 2s) 



24(1 + s)f+2 
90(1 + s)P+'^ 



{{2,p - 1)8" + (12p - 4)s3 + I2ps2) 
((3p - 5)sS + (18p - 30)s^ 



+ (45p - 75)s^ + {60p - 90)s^ + (45p - 45)5^) . 

(16) 

The variance of F, which is equal to varji?}, can be de- 
rived from this expression. Up to order of N^, we obtain 

1 s^ + 6s^ + 15s^ + 15s^ 
var{r} = — TT-r^ • (17) 



15 



(l + s)6 



This results in 



vax{g} = 



15 



(18) 



irrespective of ni in the (lifFiisivc regime of iY ^ s ^ 1. 
This indicates that the perfectly conducting channels do 
not affect var{.g} in the large- A/' limit. 

3. Monte Carlo Approach 

To obtain var{gf} without the restriction of TV 3> m, 
we employ a Monte Carlo approach based on an analytic 
expression of the probability distribution for transmis- 
sion eigenvalues. Let us introduce a set of variables Xa, 
related to Xa by Xa = sinh^ Xa- We analytically obtain a 
simple approximate expression of the probability distri- 
bution P{{xa}; s) from the exact solution of the DMPK 
equation. The DMPK equation has been solved exactly 
for the ordinary case of m = 0.^^ The exact solution 
for an arbitrary m has been obtained in ref. 20. In the 



notation which is convenient for our purpose, the exact 
probability distribution is given by 

N 

P{{xa}; s) = const. (sinh^ Xa — sinh^ Xb) 

a,5=l(a>6) 

N 



sinh 2xa sinh^ 



a=l 



X det {/6(-sinh^a;a)}^ 2,., 



,N ' 



(19) 



where det{Aab}a,b=i,2,...,N denotes the determinant of 
the N X N matrix A and 



/6(-sinh2a;„) = / dkk^^^-'^^cl,ik)e-T^' 
Jo 

X F„(A;, -sinh^ Xa) 



with 

d(fc) 



1 |r(2i±i±i^) 



47r r2(m-M)|r(ifc)|2' 

Fm{k, — sinh^ x) 

^'m + l-ik m + l + ik , . o 
= t I , , m -|- 1; — smh x 



(20) 
(21) 



(22) 

When s/AN <C 1, the dominant contribution to the in- 
tegration over k comes from the large-A: region of A; > 
^/4N/s 1. In this region, we can approximate that 

1 



22™+i(to!)2 
Fm{k, — sinh^ x) 

l-ik 1+ik 



(23) 



,m-|-l;-sinh^a; ) . (24) 



We express the asymptotic form of Fjn{k, — sinh x) in 
terms of a Bessel function. We start with the following 
expression 

^/1-ifcl-hifc , .,2 
F I — ,m-|- 1; -smh x 

/coshx\ p-m. ij2a;). (25) 
\ smh X ) 

We replace the Legendre function by its asymptotic form 



Pi";™ (cosh 2x) 



X cos I kx 



(26) 



. 7rsinh2a;^ 

ran TT^ 
~2~ ~ 4. 

for large k. This expression can be derived from 
eq. (8.723) of ref. 21. Combining eqs. (24)-(26) and 
the asymptotic form of the Bessel function Jmiv) ^ 
cos(y — rmr/2 — 7r/4) for large y, we arrive at 

/ cosh x\"^ f 2.x 



Fm{k, — sinh a;) ~ m! — 



k j \ sinh X 

^ Jm{k^)' 



\ sinh 2.T 



(27) 
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Substituting cqs. (23) and (27) into eq. (20) and carrying 
out the /c-integration, we obtain 

(6—1)! [ Xa cosh .To 



76(-sinh Xa) 



(- 



'2iXf, 



\ sinh 2xa 



sinh x„ 



^6-1 



Nxi 



(28) 

where L'^_i is the Lagucrre polynomial. Substituting 
eq. (28) into eq. (19) and using the following relation 



JV 



const. Y[ i^l - ^b) 

a,6=l(a>6) 



(29) 



wc obtain 

P{{Xahs) 



N 



const. JJ^ { (sinh^ Xa — sinh^ x^) {x^ — } 

a,5=l(a>6) 



N 



a=l 



]j[ <^ (xa sinh 2x0)™^^ e 



(30) 



For the case of to = 0, the identical probability distri- 
bution has been obtained in ref. 19. It is convenient to 
rewrite eq. (30) as 

P{{xa}; s) = const. e-^«^»>\ (31) 

where 

N 



/ 

a=l ^ 



TO + - ) In sinh 2xa 



N 

— ^In I sinh^ Xa — sinh^ xi,\ + In | — a;^ | ^ 

a,6=l(a>6) 

(32) 

with 7 = N/s. From eq. (31), we find that the average 
of VP = {J2a=i 1/ cosh^ Xa)P is expressed as 

/■OO 

{TP) = Z-^ / da;i • • • da;jvrf e-^^^^"" (33) 
Jo 



with 



POO 

Z= da^i-.-dajAre-^^^^o". 
Jo 



(34) 



Strictly speaking, our approach is justified only when 
is sufficiently small, since we have assumed s/AN <C 
1 in the derivation of eq. (30). However, we can expect 
that it is qualitatively reliable even when 7"^ is not 
small. To see this, let us consider the large-7~^ limit. 
In this limit, the variables {xa} become much greater 
than unity and are widely separated with each other, so 
that we can assume 1 <^ xi <^ X2 ^ ••• xn- Under 
this assumption, H{{xa}) can be approximated as 

JV 

H{{xa}) = J2 (^^a - (2a - 1 + vm)xa) (35) 

a=l 



with rj = 2. Equation (35) is identical to the correct 
expression in the large-7~^ limit^^ if we substitute r] = 
1. This implies that the probability distribution given 
in cq. (31) is qualitatively reliable even for not small 
7~^ although it overestimates the influence of perfectly 
conducting channels. 

Note that we can interpret H{{xa}) as the Hamil- 
tonian function of N classical particles in one dimen- 
sion. This analogy allows us to adapt a Monte Carlo ap- 
proach to numerical calculations of (F^). Using a simple 
Metropolis algorithm,^^' ^'^ we compute yax{g} for TV = 6 
and 20 with m = 0, 1, 2, 3 as a function of 7"^ = s/N. 
The results are shown in Fig. 1, where 7"^ is restricted 
tol>7~^>0 and the average for each data point is 
taken over 5 x 10^ Monte Carlo steps. Figure 1 shows that 




Fig. 1. The variance of g for (a) N = 6 and (b) A'' = 20 as a 
function of 7^^ = s/N. Circles, squares, triangles and diamonds 
correspond to m = 0, 1, 2 and 3, respectively. The solid lines are 
guide to eyes. 



vai{g} depends on to. In the ordinary case of to = 0, we 
observe from this figure that varjg} increases with in- 
creasing 7"^^ for > 7"^^ > and approximately 
takes a constant nearly equal to the universal value 1/15 
for 7~^ > N~^. However, deviation from this universal 
behavior arises when to 7^ 0. For to, 7^ 0, we observe that 
ya.r{g} does not take a constant value, but decreases with 
increasing 7"^ for 7"^ > N~^. The decrease of ya.r{g} 
becomes more pronounced with increasing to. We also 
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observe that the deviation in the case of TV = 6 is more 
noteworthy than that in the case of iV = 20. We conclude 
that the universahty of conductance fluctuations breaks 
down for m ^ unless N is very large. This conclusion is 
not inconsistent with the result of §2, because the 
expansion approach is justified only in the large- A?' limit. 

4. Summary 

We have studied the variance of the dimensionless con- 
ductance g in disordered wires with unitary symmetry 
in the diffusive regime. We have focused on the case in 
which the number of left-moving conducting channels is 
N + m, while that of right-moving ones is N. In this 
case, m left-moving channels become perfectly conduct- 
ing. First, we have analytically obtained the variance as 
a function of s = L/l hy using an N~'^ expansion ap- 
proach. We have shown that the variance is given by 
va,i{g} = 1/15 irrespective of m when N ^ maxj™, 1}. 
This indicates that we cannot detect influences of the 
perfectly conducting channels in the large- iV limit. Sec- 
ond, we have numerically calculated varjg} by using 
a classical Monte Carlo approach to study the case in 
which the ratio m/N is not very small. We have shown 
that var{(/} approximately takes a constant value nearly 
equal to 1/15 only in the ordinary case of m = 0, but 
deviation from it appears when m > 1. We have also 
shown that the deviation becomes more pronounced with 
increasing m, particularly for small N. We conclude that 
the universality of conductance fluctuations breaks down 
for m 7^ unless N is very large. 
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Appendix: Derivation of (F^) 

We substitute eqs.(ll)-(13) into eq. (6) and equate the 
coefficient of Np+i (g = 1, 0, -1) in the left-hand side 
with that in the right-hand side. We obtain 

/p,o(s) = -p/p-M,o(s), (A-1) 

= -mpfpfiis) - pfp+i,iis), (A-2) 

/p,2(s) = -mpfp^iis) + p{p - 1) (ip_i,o(s) - fcp_i,o(s)) 

-p/p+i,2(s). (A-3) 

In a manner similar to this, we substitute eqs. (11)-(14) 
into eqs. (8)-(10) and obtain 

i;,o(s) = 2/p+i,o(s) - (P + 3)jp+i,o(s), (A-4) 

^p,o(s) = 6ip+i,o(s) - {P + 5)kp+ifl{s) - 3lp+ifl{s), 

(A-5) 

lp,oi-^) = 4jp+i,o(s) -ip + 6)lp+uo{s). (A-6) 
We obtain (r^) by solving eq. (A-l)-(A-6). 



Prom eq. (A-1), we obtain 



Substituting this into eq. (A-2), we obtain 

mp 



(A-7) 



(A-8) 



If we assume that /p,i(s) = PXpi^) with Xp+iC-s) = (1 -|- 
s)~^Xp{s), the above equation is reduced to 



, , , p + l , , m 

Xp(«) + YT^Xp(^) = -(Y^- 



(A-9) 



Solving this equation with /p,i(0) = PXp(O) = 0) we ob- 
tain Xp{s) = -(m/2)(s2 + 2s)/(l -I- s)P+^. This results 
in 

We need jp_i,o(s) and kp-ifi{s) to obtain /p,2(s)- Sub- 
stituting eq. (A-7) into eq. (A-4), we obtain 

j'pfiis) + (P + 3)ip+i,o(s) = (YTl)i^- ''^'^^^ 

If we assume that jpfl{s) satisfies jp+ifl{s) = (1 -|- 
s)~^jp,o{s), eq. (A-11) is reduced to 

The solution of eq. (A-12) with jp,o(0) = 1 is given by 

We here treat lpfi{s) which is necessary to obtain kp Q{s). 
Substituting eq. (A-13) into eq. (A-6) and assuming that 
lp,o{s) satisfies Zp+i,o(s) = (1 + s)~^lpfi{s), we obtain 

^p,o(s) + Y^'P'0(*) 

= 3(1 + sy+^ + + 6. + 3) . (A-14) 
The solution of eq. (A-14) with Zp,o(0) = 1 is given by 

.o(^) = 9(i+\)p+6 (4^' + 24s^ + 60s* + 845^ 

-f 72s2 + 36s + 9). (A-15) 

We substitute eqs. (A-13) and (A-15) into eq. (A-5) and 
assume that kpfi{s) satisfies kp+ifl{s) = {l + s)~^kpfi{s). 
We then obtain the following differential equation 



: (85*^-1-488^ -I- 1205^^-1-16253 



3(1 + s)P+7 

-|-126s^-|-54s-F9). (A-16) 
The solution of eq. (A-16) with kpfl{0) = 1 is given by 
1 



kpfiis) 



15(1 + s)P+6 



Ss'^ + 485^ + 1205* -I- 1655^ 
-h 1355^-1- 605 -1- 15). (A-17) 
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Wc now turn to the evaluation of ,fp.2{s). Substituting 
eqs. (A-10), (A-13) and (A-17) into eq. (A-3), we obtain 



f'pfiis) +Pfp+l,2{s) 



9 9 

m p 



+ 



2(1 + .s)P+i 



+ 2s) 



15(1 + s)P+5 



(2s^ + 125^ + 30s* + 40s^ 

+ 30s^ + 15s). (A-18) 

If we decompose fp,2{s) as fp,2{s) = m?Up{s) +Vp{s), we 
obtain the differential equations for Up{s) and Vp{s) as 



Up(s) +pup+i{s) 



2(1 + ,s) 



(A-19) 



<(s) +V0p^,{s) = ^0-^^{2s^ + 12/ + 30s* 

+ 40s^ + 30s2 + 15s). (A-20) 

We can solve eq. (A-19) by assuming Up{s) = p^Hpi^s) + 
pvp{s) with Mp+i(s) = (1 + s)~Vp(s) and fp+i(s) = 
(1 + s)~^Vp{s). Equation (A-20) can also be solved by 
assuming that Vp{s) = p^^p{s) + pTp{s) with ^j,_|_i(s) = 
(1 + s)~^^p{s) and Tp+i(s) = (1 + s)~^Tp{s). Combining 
the resulting Up{s) and Vp{s), we obtain 



/p,2(s) 



rn?p 



24(1 + s)P+2 



+ 



((3p - 1)5* + (12p - A)s^ + 12ps^ 
■((3}5-5)s*^ + (18}j-30)s^ 



90(1 + s)P+* 

+ (45j9 - 75)s* + (60p - 90)s^ 

+ (45p-45)s2). (A-21) 



Substituting eqs. (A-7), (A-10) and (A-21) into eq. (11), 
we finally obtain eq. (16). 
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